Interplanetary and interstellar plasma turbulence 



A A Schekochihin 12 , S C Cowley 2 3 and W Dorland 4 

1 DAMTP. University of Cambridge, Cambridge CB3 OWA, UK 

2 Department of Physics, Imperial College London, London SW7 2BW, UK 

3 Department of Physics and Astronomy, UCLA, Los Angeles, CA 90095-1547, USA 

4 Department of Physics, University of Maryland, College Park, MD 20742-3511, USA 

E-mail: as629@damtp.cam.ac.uk 
Abstract. 

Theoretical approaches to low-frequency magnetized turbulence in collisionless and 
weakly collisional astrophysical plasmas are reviewed. The proper starting point 
for an analytical description of these plasmas is kinetic theory, not fluid equations. 
The anisotropy of the turbulence is used to systematically derive a series of reduced 
analytical models. Above the ion gyroscale, it is shown rigourously that the Alfven 
waves decouple from the electron-density and magnetic-field-strength fluctuations and 
satisfy the Reduced MHD equations. The density and field-strength fluctuations (slow 
waves and the entropy mode in the fluid limit), determined kinetically, are passively 
mixed by the Alfven waves. The resulting hybrid fluid-kinetic description of the low- 
frequency turbulence is valid independently of collisionality. Below the ion gyroscale, 
the turbulent cascade is partially converted into a cascade of kinetic Alfven waves, 
damped at the electron gyroscale. This cascade is described by a pair of fluid-like 
equations, which are a reduced version of the Electron MHD. The development of these 
theoretical models is motivated by observations of the turbulence in the solar wind and 
interstellar medium. In the latter case, the turbulence is spatially inhomogeneous and 
the anisotropic Alfvenic turbulence in the presence of a strong mean field may coexist 
with isotropic MHD turbulence that has no mean field. 
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1. Introduction 

Rapid progress in astronomical instrumentation has made it possible to observe 
astrophysical plasmas with ever greater spatial resolution. This has allowed astronomers 
to probe not only the bulk, large-scale motions and fields but also to measure, either 
directly or via line-of-sight integrated quantities associated with the emission and 
propagation of light, the small-scale fluctuations of plasma velocity, density, magnetic 
and electric fields. These turbulent fluctuations existing in a broad range of scales are 
a common property of astrophysical plasmas. While astrophysical turbulence occurs 
in a variety of vastly differing conditions, its physical characterization is based on a 
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number of universal features. In most cases, the source of energy is in the form of 
random stirring or instabilities associated with the scale of the astrophysical object of 
interest. The energy injected at large scales cascades to much (typically many orders 
of magnitude) smaller scales to be dissipated into heat. A signature property of the 
turbulent cascade that connects these vastly disparate scales is power-law spectra of the 
fluctuating quantities. These have been observed in the solar wind (SW), e.g. [351 12*1 12"^] . 
the interstellar medium (ISM) S a l ax Y clusters [S31EZI, etc. In all of the cited 

examples, the reported spectra had, or were consistent with, Kolmogorov scaling k~ 5 ^ 3 . 

In this paper, we shall concentrate on the SW and ISM and outline both the 
qualitative understanding that currently exists of the turbulence in these media and 
a formal mathematical description of this turbulence that must underlie the future 
analytical and numerical investigations of it. The proper starting point for such a 
description is the kinetic plasma theory because the turbulent plasmas we are interested 
in are either collisionless (in the SW, the particle mean free path is comparable to the 
distance from the Sun to the Earth) or only weakly collisional, meaning that the mean 
free path A m f p exceeds the ion gyroradius pi (in the ISM, A m f p ~ 10 12 cm, p { ~ 10 9 cm). 

In many cases, it is plausible to think of plasma turbulence at scales much smaller 
than the energy- injection scale as an ensemble of interacting MHD waves propagating 
along a dynamically strong background magnetic field (the mean field) associated with 
the large scales |H2|- Goldreich and Sridhar j20] (henceforth, GS) conjectured that in 
such a turbulence, (i) all electromagnetic perturbations are strongly anisotropic, so that 
the characteristic wavenumbers along the field are much smaller than those across it, 
fcii <C k±; and (ii) the interactions between the Alfven waves are strong, i.e., the Alfven 
time and the nonlinear interaction time are comparable to each other: 

uj ~ k\\v A ~ k±u±, (1) 

where to is the typical frequency of perturbations, v A is the Alfven speed, and u±_ is 
the velocity fluctuation perpendicular to the mean field. This assumption, known as 
the critical balance, removed dimensional ambiguity from the MHD turbulence theory 
and led to the Kolmogorov scaling of the Alfven-wave energy spectrum, kj 5 ^ 3 and to 
the relation k\\ ~ k 2 / 3 (for a historical review, see [201; in | Appendix A[ we give a brief 
outline of the GS theory and related scaling arguments for MHD turbulence). 

The anisotropy of MHD turbulence is supported by observations of the SW the ISM 
(see reviews j2HEH]) and by numerical simulations P^IT^HHl ll2j. In what follows, this 
anisotropy emerges as the key simplifying feature used to derive a reduced version of 
the plasma kinetic theory that describes low-frequency MHD turbulence. This is done 
in ^21 where our exposition is motivated by the observations of the collisionless SW. We 
show how the descriptions known as Reduced MHD, Kinetic MHD, Electron MHD and 
Gyrokinetics fit into a single theoretical framework. In JJJ we explain how the same 
approach works for the turbulence in parts of the ISM and how this type of turbulence 
differs from the isotropic MHD turbulence, which does not have a mean field. We argue 
that the latter kind of turbulence may also be present in the ISM and in galaxy clusters. 




Figure 1. Spectra of electric and magnetic fluctuations in the SW — adapted with 
permission from figure 3 of [2] (copyright 2005 by the American Physical Society). 
The bold dashed (red) lines are reference slopes added by us. We also inserted 
"KMHD" , "GK ions" and "ERMHD" to indicate the scale intervals where these 
analytical descriptions are valid: k± <C p^ 1 (sec H2.3|) . k± ~ p^ 1 (see ^2.6[) and 
p^ 1 <C k± <C Pe 1 (sec H2.70 . respectively. 



2. Solar wind and the collisionless MHD turbulence 

Spectra of electromagnetic fluctuations in the SW extend across a broad range of 
collisionless scales. Above the ion gyroscale (pi ~ 100 km), the spectra of the electric and 
magnetic field measured by spacecraft at 1 AU from the Sun fit the k~ 5 ^ 3 law and follow 
each other with remarkable precision (figure Q). Since the electric field is directly related 
to the plasma velocity (at scales above pi, it is the E x B drift velocity), this can be 
interpreted as a signature of Alfvenic turbulence. How do we describe such a turbulence 
at collisionless scales? Let us first consider scales larger than the ion gyroscale. 

Note that the plasma beta $ is taken to be order unity in what follows, as is 
appropriate both for the SW and the ISM. It is useful to remember that in this regime, 
the ion inertial scale is comparable to pi. 



2.1. Kinetic MHD 



For kpi <C 1, the magnetic field impedes free particle motion across the field lines and 
the kinetic theory reduces to the so-called Kinetic MHD (KMHD) [HI], which has most 
features of the MHD description, but allows for anisotropic pressure: 
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— = B Vit - BV ■ u, (4) 

dt v 1 

where p is mass density, u velocity, B magnetic field, b = B/B, and d/dt = d/dt+u-'V. 
The pressure tensor is calculated kinetically: p± = J2 s m s f d i v{y 2 _/2)f s and p\\ = 
^2 s m s J d 3 v(v\\ —u\\) 2 f s , where the distribution function f s (t,r,v±,v») satisfies 

Df s , lDBv x df s ( t Du ± , v\b-VB q,E\\\df a . _ 

0, (5) 



B Dt 2 dv± \ Dt 2 B m s J dv u 

where D / Dt = d/dt + (u±_ + vn b) ■ V and u± = u — u- bb. In the above, m s and q s = ±e 
are the mass and charge of the particles of species s (ions, electrons). The parallel 
electric field E» is determined from the quasineutrality condition Y2 S Qs^s = 0, where 
n s = J d 3 v f s (the number density). Note that p = miUi and un = / d 3 vv\\fi, so 

(0) and the parallel component of (j3J) can be derived from (J3J). 

We consider a uniform static equilibrium with a straight mean field in the z 
direction, so B = B z + 5B and p = po + Sp, p± = po + 5p±, p\\ = po + Spn, where B , 
Po, Po are constant in space and time. 



2.2. The ordering 

The anisotropy of the turbulence allows us to systematically expand ©-© in e ~ 
k\\/k_\_. The key step in setting up such an expansion is to estimate the strength of 
the fluctuations by adopting the critical-balance conjecture (0), but as an ordering 
assumption rather than a detailed scaling prescription: this means that the wave 
propagation terms are assumed to be same order as the nonlinear interaction terms 
(the turbulence is strong). This leads to the following ordering: 

5p u± 5p± 5p\\ 5B± SB\\ 

Po v A v A po Po Bo Bo 
where va = -Bo/V^Po- Two auxiliary ordering assumptions have been made: (i) the 
perpendicular velocity and magnetic-field fluctuations are Alfvenic (5B±/B ~ u±/va)', 
(ii) the Alfvenic fluctuations are same order as magnetic-field-strength (SB\\), density 
and pressure fluctuations (in the collisional MHD limit, these correspond to the slow 
waves and the entropy mode). The validity of the latter assumption depends on how 
the turbulence is stirred. In astrophysical contexts, the large-scale energy input may be 
assumed to inject comparable power into all types of fluctuations. We also assume that 
the characteristic frequency of the fluctuations is u ~ knVA- The fast waves are, thus, 
ordered out, because their frequency is ~ fcyiu/e. 



2.3. Alfven-wave turbulence: Reduced MHD 

The Alfvenic fluctuations are two-dimensionally solenoidal: since V ■ u = 0(e 2 ) [from 
(J2J)] and V • SB = exactly, separating the 0(e) part of these divergences gives 
Vj_ • u± = • 5B±_ = 0. Therefore, to lowest order, we may write u± = z x and 
5B±/ y/Anpo = z x Vi$. Equations for the scalar fields $ and \& (the stream and flux 
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functions) are obtained by substituting these expressions into the perpendicular parts of 
(J3J) and (jlj) — of the former the curl is taken to annihilate the pressure term. Keeping 
only the terms of the lowest order, 0(e 2 ), we get 

^vi$ + {$,vi$} =v A ^vl* + {*,V i ji}, (7) 

d d 

+ = WA _$, (8) 

where {$, tf/} = z ■ (Vi$ x Vj_\l/) and we have taken into account that, to lowest order, 

d 9 n 9 i L ^ 9 SB 1- « 9 1 r T 

- = - + ^.V, = - + { $,..},6.V = - + ^-.V, = - + -R...}.(9) 

The closed system known as the Reduced MHD (RMHD), was derived originally 

from the fluid MHD equations for the studies of stability of fusion plasmas (HElEni- We 
have now shown that the Alfven-wave cascade in a collisionless plasma is described by 
the RMHD equations (0)-(Hj) oil the way down to the ion gyroscale. The Alfven waves 
are decoupled from the other fluctuation modes: density and magnetic-field-strength 
fluctuations (slow waves and entropy fluctuations in the collisional limit; cf. fZT\ 13 7j). 
Introducing Els asser fields £ = $ ± we may rewrite (f7|) -([B|) as follows 

| ViC± T VA ^ VlC± = ~\ ^ C+ ' ViC ^ + ^ ViC+ ^ T Vl ^ C+ ' C ^ • (10) 
Thus, the RMHD, like the MHD, supports wave packets of arbitrary shape and 
magnitude propagating in one direction at the Alfven speed va- if C = or C + = 0, 
the nonlinear terms vanish and the exact solution for the other Elsasser potential is 
= f ± (x, y,z=f vxt), where is an arbitrary function. The Alfven-wave cascade is 
a result of interactions between counterpropagating wave packets [3~2j . 

It is this Alfvenic component of the plasma turbulence to which the GS scaling 
theory of MHD turbulence flAppcndix A| ) applies. In the SW, it is observed via 
in situ measurements of the fluctuating magnetic and electric fields j2] (see figure 
HJ. The latter directly probe the velocity fluctuations because, to lowest order in e, 
u± = cExB/B 2 = (c/Bo)zx V_l</>, where <p is the scalar potential. Clearly, $ = c(f)/B . 



2.4- Density and magnetic-field-strength fluctuations: Kinetic Reduced MHD 

In order to determine 5n e and SB», we must use the kinetic equation (0). 2 
The lowest-order (equilibrium) distribution is taken to be a Maxwellian: F 0s = 
noe~ v2 ^ v ^=/(7iv 2 hs ) 3 ^ 2 , where v t hs = (2T /m s ) 1,/2 is the thermal speed of species s and T 
is temperature. 3 We let f s = F 0s + 5f s , where Sf s /F 0s ~ e and apply the ordering © 
to the kinetic equation (jSJ). 

2 These quantities cannot be derived from J5J) and the parallel part of (0} because (i) 10 has already 
been used to determine V ■ u, a 0(e 2 ) quantity; (ii) the parallel part of contains , whose evolution 
equation, the parallel part of requires Sp± — Spn, so it|| can only be calculated kinetically. 

3 The assumption of an isotropic equilibrium was implicit when we adopted an isotropic zeroth-order 
pressure po at the end of H2.ll Strictly speaking, in a collisionless plasma such as the SW, the equilibrium 
distribution does not have to be Maxwellian or isotropic. The conservation of the first adiabatic 
invariant, ji = v\j1B, suggests that temperature anisotropy with respect to the magnetic-field direction 
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The electron kinetic equation can be further simplified by a subsidiary expansion 
in (m e /m i ) 1 '' 2 j^HJ. To lowest order, 

b-V6f e + ^E l{ F 0e j = 0. (11) 

Since f d 3 v Sf e = 5n e , the inhomogeneous solution of this equation is Sf e = (5n e /no)F 0e 
(the electrons are isothermal). The homogeneous solution satisfies b ■ V5/ e = 0, i.e., it 
is constant along the perturbed field lines and is constant everywhere if the field lines 
are assumed to be stochastic. Thus, E\\ = — (T /eno)b ■ V5ri e . Substituting this into 
the ion kinetic equation, we have, to lowest order, 0(e 2 ), 

d / „ „ v 2 , SB\\ „ \ ? _ / „ „ 5n 



T A5fi- -^-^F 0i +v l{ b-V 5ft + — F 0l = 0. (12) 

Finally, we calculate 5n e and 5B\\. From quasineutrality, Sn e = Srii, so 
Sn e 1 



d'vSU (13) 

n n J 

To calculate SB», we first revisit the the perpendicular part of 0- In the lowest order, 
0(e), it reduces to the perpendicular pressure balance: Vj_ (5p± + B 5B^/An) = 0, 
whence 5B\\ = —(47t/Bq)5p± (this is why the fast waves disappear under our ordering). 
Now 5p± = 5p± e + 5p±i. Using 5f e = (5n e /no)F 0e to get Spj_ e = T 5n e , equation (fT3j) to 
express Sn e , and calculating 5p±i from 5 fa, we find 



,2 



^ = 4±/Wl + ^U/„ (14) 

Bo 2 n J V <i/ 
where $ = StihoTo/Bq = v 2 hi /v\. Note that u\\ = (1/tio) / d 3 vv\\5fi is not required to 
solve the equations, but can be calculated from the solution. 

Together with (|T^ |l -(|Ti |l form a closed system that describes the anisotropic 

turbulence above the ion gyroscale in a collisionless magnetized plasma. We shall refer 
to this hybrid fluid-kinetic theory as Kinetic RMHD (KRMHD). The nonlinearity enters 
in (|T2*|) via the derivatives defined in © and is due solely to interactions with Alfven 
waves. Thus, the cascades of density and magnetic- field- strength fluctuations occur via 
passive mixing by Alfven waves, with no energy exchange (cf. |2*H I57|). 



2.5. Parallel and perpendicular cascades 

Let us transform ()12|) to the Lagrangian frame associated with the velocity field u± 
of the Alfven waves: (t,r) — > (t,r ), where r(t,r ) = r + f Q dt' u±(t', r(t', r )). In 

(Tj_o 7^ Tj| ) may exist. Such anisotropy gives rise to several high-frequency plasma instabilities 
and it is plausible to assume that fluctuations associated with them will scatter particles and limit the 
anisotropy (e.g., [HO])- While there is no definitive analytical theory quantifying this idea, it has some 
support in the SW observations that indicate that the core particle distribution is only moderately 
anisotropic 0Oj- We believe, therefore, that assuming a Maxwellian equilibrium is an acceptable 
simplification. We also take Xbi = Tq c (generalising to ^ Tq c is straightforward). Note that in 
plasmas such as the ISM, where collisions are weak but non-negligible f H3.1J) . the Maxwellian equilibrium 
is rigourously justifiable if the ion collision rate is ordered vu ~ u> within the e expansion |26| . 
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this frame, d/dt [defined in (0)] becomes d/dt. Equation has the Cauchy solution: 
B(t) = [p(t)/p{0))B(0) ■ V r, where V = d/dr . Then b ■ V = 6(0) • (V r) ■ V = 
6(0) ■ Vo = d/dlo, where Iq is the arc length along the magnetic field line taken at 
t = [if 5B±(0) = 0, l = zq\. Thus, in the Lagrangian frame associated with the 
Alfven waves, (fT^jl is linear. It does not, therefore, support a cascade of 5n e and 5B\\ to 
smaller scales parallel to the perturbed magnetic field, i.e., 6 ■ V of these fluctuations 
does not change with time. In contrast, passive mixing by the Alfven waves does cause 
a perpendicular cascade of 5n e and 5B\\ — i.e., a cascade in k±. 

Unlike (|T2|) . the RMHD equations ©-© in the Lagrangian frame do not reduce 
to a linear form, so the Alfven waves should develop small scales both across and along 
the perturbed magnetic field. The scale-by-scale critical balance (0) conjectured by GS 
leads to the relation k\\ ~ k 2 / 3 (see | Appendix A ) . 

Using the linearity of ()12j) in the Largangian frame, it is straightforward to show 
that density and field-strength fluctuations are damped. The dispersion relation is 

]k^z z (sfc) = ~ 2 i 1 ~ ^ ± ' (15) 

where Z is the plasma dispersion function and c^o and fcyo are the Lagrangian frequency 
and wave number (k\\ ~ 6 • V). When fa ~ 1, all solutions of (fTo]) have damping rates 
Im(u;o) ~ — |fc||o|^thi ~ — |^||o| w a- 4 If no parallel cascade of 5n e and SB\\ develops, the 
parallel wavenumber k\\ of these fluctuations with a given k± does not grow with k±, 

2 /3 

so, for large enough k±, it is much smaller than the parallel wave number k\\A ~ kj_ of 
the Alfven waves that have the same k±. This means that the damping rate is small 
compared to the characteristic rate k\\A v A a t which the Alfven waves cause 5n e and 8B\\ 
to cascade to higher k±. One is then led to conclude that, despite the kinetic damping, 
5n e and SBu should have perpendicular cascades extending to the ion gyroscale. 

The validity of this conclusion is not quite as obvious as it might appear. Lithwick 
and Goldreich [37] argued that the dissipation of 5n e and 8B\\ at the ion gyroscale 
would lead these fluctuations to become uncorrelated at the same parallel scales as 
the Alfvenic fluctuations by which they are mixed, i.e., k\\ ~ k\\A- The damping rate 
then becomes comparable to the cascade rate, causing the cascades of density and field- 
strength fluctuations to be cut off at fc[|A m f p ~ 1. In the SW, this would mean that no 
such fluctuations should be detected above the ion gyroscale. Observational evidence 
is at odds with this conclusion: the density fluctuations appear to follow a k~ 5 ^ 3 law 
at kj_pi 1 [TT] . as they should if they are passively mixed and not damped (see 
Appendix A| ) . The same is true for the fluctuations of the field strength jHl • It is not 
clear why Lithwick and Goldreich's argument fails, but it is, perhaps, useful to point 
out two potential pitfalls: (i) in order for the dissipation terms, not present in ([12 )1 - (|14 |) . 

4 For (3i ^> 1, the weakest-damped solution is uio ~ — «|fc|| |i>A/-\/ 7r /^- This is the anisotropic limit 
(ku/k± <C 1) of the more general effect known as Barnes, or transit-time, damping Note that we 
carried out the expansion in small fc||/fc_i_ before taking the high-/? limit. A more standard approach in 
the linear theory of plasma waves is to leave k\\/k± arbitrary and take the high-/3 limit first 
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to act, the density and field-strength fluctuations should reach the ion gyroscale in the 
first place; (ii) the damping rate of these fluctuations, even if km ~ fciu, is never much 
larger than the cascade rate, so it may be necessary to have a quantitative calculation 
of the interplay between the kinetic damping, mixing and the dissipation at k±pi ~ 1 
in order to determine the efficiency of the cascade. 

2.6. Gyrokinetics 

At k±pi ~ 1, the KMHD description breaks down and the Alfvenic fluctuations are no 
longer decoupled from the kinetic component of the turbulence. They are mixed with the 
fluctuations of the density and magnetic-field strength and dissipated via the collisionless 
damping discussed in fc |2.5l — the observed flattening of the density-fluctuation spectrum 
as k±pi approaches unity is likely to be due to this energy exchange with the Alfven 
waves. The damping leads to ion heating, an astrophysically interesting problem in its 
own right, e.g., in the theories of coronal heating ^3] and accretion discs [IH]. The 
amount of heating suffered by the ions is a nontrivial issue because only part of the 
turbulent energy is dissipated at kj_pi ~ 1. The rest is converted into a cascade of kinetic 
Alfven waves (KAW) that extends to the electron gyroscale — a feature observed in 
the SW j2El E] . Quantitative theory or numerical modeling of the energy dissipation 
and conversion processes at k±pi ~ 1 can only be done in the fully kinetic framework. 
However, the anisotropy of the fluctuations leads to a substantial simplification of the 
full plasma kinetic theory. If the ordering based on the assumptions of anisotropy and 
critical balance ( A2.2JI is applied, the plasma kinetics reduce to gyrokinetics (GK) — a 
low-frequency limit well known in fusion science jl()j . A simple derivation of GK based 
on the ordering of £ 12.21 is given in along with a detailed GK treatment of the 
linear collisionless damping at k^pi ~ 1. The GK is a valid approximation at all scales 
that are of interest in the context of low-frequency astrophysical turbulence, down to the 
electron gyroscale and below. This broad range of validity and the long experience of GK 
simulations developed in fusion research make GK an ideal tool for numerical modeling 
of astrophysical turbulence 5 and a good starting point for analytical theory. 

The RMHD and KRMHD equations (S 92.3I2.4|) can be derived from GK by means 
of two subsidiary expansions: first in (m e /mj)^ 2 , then in k±pi <C l. 6 This and various 
other limits of the GK description of turbulence in weakly collisional astrophysical 
plasmas are worked out in [H2]- While the k±pi ~ 1 regime requires solving the ion GK 
equation (electrons remain isothermal the KAW turbulence at p^ 1 <C k± <C p" 1 

is described by another well known fluid model, the Electron MHD (EMHD) 

5 A programme of such numerical studies, using the GS2 code http: / /gs2.sourceforge.net7| , is currently 
underway (supported by the US DOE Center for Multiscale Plasma Dynamics). 

6 This means that the k\\/k±_, (rrie/mi) 1 / 2 and k±pi expansions commute: KRMHD can be arrived at by 
either of the two routes: full kinetics — > k^pi expansion — > KMHD — > ku/k± expansion — > (m^/riii) 1 / 2 
expansion — * isothermal electrons — > KRMHD [this paper] or full kinetics — > k^/k± expansion — > GK 
[221 — > (me/mi) 1 / 2 expansion — > isothermal electrons — > k±pi expansion — > KRMHD |52j . 
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2. 7. Kinetic Alfven waves: Electron Reduced MHD 

When k±pi 3> 1, the ions are unmagnetized and have a Boltzmann distribution: 
fi = Foi(v) exp (— e0/T o ), where is the scalar potential. The electrons are magnetized 
(k±p e <C 1) and can be shown to be isothermal in essentially the same way as in 
where fill)) is still valid to lowest order in (m e /mj) 1//2 . Then 

Sn e _ 5rk _ _e0 _ 

n n T A .Bq ' 

The last equality follows from the (perpendicular) pressure balance similarly to the way 
it was done in < j2.4l using 5p±i = 5p± e = Sn e To (the ordering of < J2.21 which eliminates 
the fast waves, continues to be valid). The EMHD equations now follow from the density 
and parallel velocity moments of the electron kinetic equation, which is similar in form to 
f)12|) . The rigourous GK derivation is given in [52]. Here, we adopt a more conventional 
approach by noting that in the limit kj_p e <C 1, the magnetic field is frozen into the 
electron fluid and satisfies (@J) with u replaced by the electron flow velocity u e |31j : 
OB 

— = V x (u e x B) = -u e ■ VB + B ■ Vu e - BV ■ u e . (17) 

We set 5B/B = (1/va)z x Vj_\1/ + z5B\\/Bq and expand (fTTj) using the ordering of 
£ 12.21 To lowest order in the k±pi ^> 1 expansion, u e can be found by taking the ions to 
be immobile and using Ampere's law: u e = Ui — j/eno = — (c/47reno) Vj_ x 5B. In the 
last term in (|17)) . the next-order compressible part of u e is calculated via the electron 
continuity equation: V-i^ = —(d/dt + u e -'V±)Sn e /n . Finally, using (fTBj) and denoting 
$ = c<p/Bo = {cTq/ 'eBo/3i)5B\\/ Bq [see (jlfij) ]. we find that the parallel and perpendicular 
components of (|17j) take the following form 

— = 2^b-V$, (19) 

where b ■ V is defined in (jUJ). We shall refer to this system as Electron Reduced MHD 
(ERMHD) — the anisotropic limit of EMHD. 7 

ERMHD describes the cascade of kinetic Alfven waves (KAW), whose linear 
dispersion relation is u> = ±k\\VAk±pi/y/T + Pi with eigenf unctions $=Ffcj_pj^//2\/l + 
To understand the nonlinear cascade, one may follow the spirit of GS theory, assuming 
anisotropy (fcy k±) and strong interactions |3E]. This argument, reviewed at the end 
of Appendix A[ leads to a kj 7 ^ 3 spectrum of magnetic fluctuations (note that for KAW- 



1 /3 

like fluctuations, SByB ~ $/ ' PiVa ~ k±^/vA ~ 5B±/B ) and to the relation k\\ ~ k[ , 
quantifying the anisotropy. Both of these scalings have been confirmed by numerical 
simulations of EMHD [12]. Note that the electric- field fluctuations in this regime 
should have a k^ z spectrum because 5E ~ kj_<f) ~ kj_pi(vA/ / c)5B . Measurements of 

7 Equation {T7} with u e = — (c/4ireno) x SB and V • u e = (the incompressible limit valid if 
Pi ^> 1) is what is normally understood by EMHD. In (|18fl - l|19|) . Pi 1S arbitrary, i.e., the electron fluid 
is not assumed to be exactly incompressible. 
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the spectra of SB and SE in the SW appear to corroborate these arguments [2] (see 
figure HJ) . 

The anisotropic KAW cascade is terminated at k±p e ~ 1 by the electron collisionless 
damping. The proper description of this process is again gyrokinetic 



3. Interstellar medium and the two regimes of MHD turbulence 

3.1. Weakly collisional limit 

The anisotropic MHD turbulence in extrasolar plasmas is largely similar to the 
turbulence in the SW. The best studied of these plasmas is the interstellar medium 
(ISM), a hot low-density plasma (no ~ 1 cm -3 , To ~ 10 4 K for the Warm ISM phase) 
that makes up most of our and other galaxies' diffuse luminous matter. Turbulence in 
the ISM is stirred by colliding shock waves caused by supernova explosions, with the 
estimated injection scale L ~ 100 pc ~ 10 20 cm jH]. One important difference with the 
SW is that in the ISM, A m f p ~ 10 12 cm is substantially smaller than L, although it is still 
larger than pi ~ 10 9 cm. Thus, collisions have to be allowed for. This can be done by 
keeping a collision integral in the GK equations and ordering the ion-ion collision rate to 
be comparable to the fluctuation frequency, ua ~ uj [2HlE2j- Both the RMHD equations 
(HJ)-© above the ion gyroscale and the ERMHD equations ()18 |) -fll9 |) below it can then 
still be derived rigourously [H2] . Collisions do not appear in these equations: in (JZj), this 
is because in the k±pi <C 1 limit the collisional transport is parallel to the field lines j^j; 
in ©, the collision terms, which give rise to Ohmic resistivity, are ordered out via the 
subsidiary expansion in [m e /mi) 1 / 2 ] the latter is also true for the ERMHD equations 
(HHJl-dini)- The kinetic part of the KRMHD system, (fT2 ]) -([n |) . remains intact except 
that the ion-ion collision integral appears in (j!2)l Thus modified, the KRMHD 

constitutes a description of anisotropic plasma turbulence valid both in the collisional 
and collisionless regime. 8 To lowest order in k±pi, the collision integral has no spatial 
derivatives, so ()12)) is still linear in the Lagrangian frame of the Alfven waves and the 
discussion of the cascades of Sn e and 5B\\ given in £ 12. 51 continues to apply. The only 
difference is that there is also collisional damping of these fluctuations, which, like the 
collisionless damping, depends solely on the variation of 5n e and <5-B|| along the perturbed 
magnetic-field lines. Indeed, in the collisional limit fc||oA m f p <C 1, (fT^ - ffTlj) reduce to a 
set of fluid equations via the standard Chapman-Enskog expansion procedure 



d 5B\\ ; d 5n e , . 

— -r 1 = b ■ Vm + -, 20 

dt B 11 dt n 

d 0? _SB 



Mil = V 



dt- ^-V^ + ^.V^.V^j, (21) 

8 Strictly speaking, this is only true for fe||A m f p (me/mi) 1 / 2 . At longer parallel scales, the electrons 
are adiabatic, rather than isothermal, ST e = STi, and the standard fluid MHD theory applies. With 
the ordering of the equations for the passive part of the turbulence are the same as (I2UH - H22II . 
except now Sn e /n = —5Ti/T — (l/(3i)SB^/B [3011221 an d the transport terms are more involved 
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d« 1 = 2dfc i h /. m\ 

dt T 3 dt n 11 V T o / 
where i/y ~ fty ~ fthi^mfp are the parallel viscosity and thermal diffusivity. The ion 
temperature is related to Sn e and 8B» via pressure balance, which is written in the form 
Sn e /n = —(l/2)ST i /T — (1/ (3i)8B\\/ B Q . Equations pU|) - P^j) describe passive cascades 
of slow waves (u\\ and SB\\) and entropy fluctuations (d~n e and SB\\) mixed by Alfven 
waves jHZj via the nonlinearities contained in d/dt and b- V and damped by anisotropic 
diffusion, which occurs purely along the perturbed magnetic-field lines. 

The observational evidence is less exhaustive for the ISM than for the SW. The 
magnetic fluctuation spectra, inferred from the structure functions of the Faraday 
rotation measure, appear to be consistent with the k ± 5 ^ 3 scaling [12], although the 
accuracy of the measurements is not high. The electron-density fluctuations, measured 
by a variety of methods, are anisotropic and also seem to have a Kolmogorov scaling 
across the entire range from L ~ 10 20 cm to pi ~ 10 9 cm — this is sometimes called 
"The Great Power Law in the Sky" P3 EH]. 9 Note that while the density- fluctuation 
spectrum appears to extend to the ion gyroscale, the scale separation between A m f p and 
Pi is not sufficient in the ISM (unlike in the SW) to distinguish this from a cutoff at 
&||A m fp ~ 1 (see fl2.5|l . which, using the GS relation k\\ ~ /c^L -1 / 3 , would imply the 
perpendicular cutoff scale ~ 10 8 cm |3*7j . 

3.2. Inhomogeneously turbulent ISM: spiral arms vs. interarm regions 

It is, in fact, simplistic to view the ISM as a homogeneous plasma. The ISM is 
a spatially inhomogeneous environment consisting of several phases (of which Warm 
ISM is one) that have different temperatures, densities and degrees of ionization ^H] 
(and, therefore, different degrees of importance the neutral particles and the associated 
ambipolar damping effects have [HI]). While the role of the molecular properties of the 
multiphase ISM is left outside the scope of this paper, we would like to discuss briefly 
another aspect of the ISM's spatial inhomogeneity: the fact that it is inhomogeneously 
turbulent. One of the most prominent spatial features of our and many other galaxies 
is the spiral arms. They are denser than the interarm regions (interarms) by a factor of 
a few [IE] and observed to support stronger turbulence [H], which is not surprising as 
the concentration of supernovae is higher. Observations of magnetic fields in external 
galaxies show that the spatially regular (mean) fields are stronger in the interarms, while 
in the arms, the stochastic fields dominate [I]. A recent study of the rotation- measure 
structure fuctions in our Galaxy [22] revealed that in the interarms, the magnetic energy 
is large-scale dominated and the structure functions are consistent with Kolmogorov-like 
negative spectral slopes, whereas in the arms, the structure functions are flat down to 
the resolution limit, meaning that the magnetic energy resides at much smaller scales 
than in the interarms. With these results in mind, let us recall that there exist two 

9 There is, however, some evidence of a k ± 3 ^ 2 spectrum as well |55| — see discussion of MHD turbulence 
scalings and the polarization- alignment theory [8] in | Appendix A| 
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asymptotic regimes of MHD turbulence, depending on the relative magnitude of the 
mean and fluctuating fields, 5B Tms /B : 

I. Anisotropic Alfvenic turbulence. This is the type of turbulence discussed so 
far in this paper. It requires that a strong mean field B is present. The turbulent 
fluctuations are much smaller than the mean field: SB rms <C Bq, u rms <C va (see 
A2.2JI . The fluctuations are Alfvenic and have a Kolmogorov spectrum, with velocity 
and magnetic fields in scale- by-scale equipartition (see Appendix A ). 

II. Isotropic MHD turbulence. In this case, no mean field is present, i.e., Bq <C 
SB rms . The dynamically strong stochastic magnetic field is a result of saturation of the 
small-scale dynamo — amplification of magnetic field due to random stretching by the 
turbulent motions. Both the small-scale dynamo and its saturation are reviewed in |5Uj . 
While the definitive theory of the saturated state remains to be discovered, both physical 
arguments and numerical evidence |HJ EH] suggest that magnetic field is organized in 
folded flux sheets/ribbons. The length of these folds is comparable to the stirring scale, 
while the scale of the field-direction reversals transverse to the fold is determined by 
the dissipation physics: in MHD with Laplacian viscosity and resistivity operators, it is 
the resistive scale. 10 The structure functions of such magnetic fields are flat (SHj, with 
magnetic energy dominantly at the reversal scale. While Alfven waves propagating 
along the folds may exist |HJ EU] , the presence of small-scale direction reversals means 
that there is no scale-by- scale equipartition between velocity and magnetic fields. 

It is tempting to explain the difference between the magnetic-field structure in the 
arms and interarms by classifying the MHD turbulence in the interarms as anisotropic 
(I) and in the arms as isotropic (II). The observational evidence cited above lends 
qualitative support to this idea and so do numerical simulations of an inhomogeneously 
turbulent MHD fluid [2E|. The turbulence in the arms should be closer to the isotropic 
variety and in the interarms to the anisotropic one for a number of conspiring reasons: 
(i) M rms is larger in the arms, so 5B Tms /B Q ~ u Tms /vA should be larger; (ii) the presence 
of the spiral mean field in galaxies is usually attributed to some form of mean-field 
dynamo [SB] and it is possible to argue plausibly that this mechanism produces stronger 
mean fields in the interarms than in the arms |54j : (iii) the mean field should be pushed 
out of the more turbulent region (arms) into the less turbulent one (interarms) by the 
diamagnetic effect of turbulence [SHI EH EE] • 

Finally, we mention another class of weakly collisional astrophysical plasmas where 
isotropic MHD turbulence is believed to exist: the intracluster medium (ICM) of the 
galaxy clusters. The turbulence in these intergalactic plasmas, which constitute the 
majority of the luminous matter in the Universe, has, in recent years, been increasingly 
accessible to observational astronomy [SHI El] • For further information, the reader is 
referred to our review 



10 In weakly collisional astrophysical plasmas, such a description is not applicable and the field reversal 
scale is most probably determined by more complicated and as yet poorly understood plasma dissipation 
processes; below this scale, an Alfvenic turbulence of the kind discussed in |2 ma y exist 
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Appendix A. Scaling theories of Alfven-wave turbulence: a brief review 

Goldreich-Sridhar turbulence 

Here we outline the key steps in the GS theory of anisotropic MHD turbulence [20J. 
A more leisurely historical review, which explains how the GS argument is related to 
the earlier (isotropic) theory of Iroshnikov j2Zj and Kraichnan j32j and to the weak- 
turbulence treatment j2TJ E] , can be found in [5U] (see also [TH] and [3B] ) • 

As in the Kolmogorov-Obukhov theory of turbulence, it is assumed that the cascade 
of energy is local in scale space and the flux of energy through scale A in the inertial 
range is scale-independent: 

— ±A ~ £ = const, (A.l) 

where e is the Kolmogorov flux, the subscript A indicates fluctuations associated with 

the perpendicular scale A, and t\ is the cascade time. It is now assumed that the 

turbulence is strong, i.e., that the Alfvenic linear propagation terms are comparable to 

the nonlinear terms: 

d v A 5u ±x , _ 

va-z u± ■ V ± <^> — , (A. 2) 

oz A 

This is the critical-balance conjecture, applied scale by scale. It is further assumed that 
the cascade time is the same as the Alfven time: t\ ~ l\\\/v A - Together with (|A.l|) - (jA.2l) . 
this immediately implies 

/ / \ 1/2 / 3 \ 1/3 

6u±x ~\t) ~ (eA)1/3 ' ^~(t) a2/3 - ( a - 3 ) 

—5 /3 

The first of these scaling relations is equivalent to a k ± spectrum of kinetic 
energy, 11 the second quantifies the anisotropy by establishing the relation between the 
perpendicular and parallel scales. The fluctuations are Alfvenic, so 5B±\ ~ 5u±_\y/ATxpQ. 

11 In terms of parallel wavenumbers. i|A.3|) means that the spectrum scales as fc^ 2 . Remarkably, recent 
SW data analysis confirms this power law |25j . 
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Polarization alignment 

While the GS theory has acquired the status of the accepted view, the failure of the 
numerical simulations [HH1 HHj to reproduce the kj 5 ^ 3 spectrum has remained a worrying 
puzzle. The numerical spectra are closer to kj 3 ^ 2 , but cannot be explained by the 
Iroshnikov-Kraichnan theory [23 E2] because the fluctuations are definitely anisotropic. 
Recently, Boldyrev [H] proposed a scaling argument that allows an anisotropic Alfvenic 
turbulence with a k ± spectrum. It is based on the conjecture that u± and 5B± align 
at small scales, an idea that has had some numerical support [SHI El E]- The alignment 
weakens nonlinear interactions and alters the scalings. 

The fluctuations are assumed to be three- dimensionally anisotropic: the three 
characteristic scales are the parallel scale l\\ along Bq, the displacement £j_ of the fluid 
element perpendicular to Bq and parallel to 5B±, and the scale A of the variation of u± 
and d~B± perpendicular both to B and to 5B^. The nonlinear terms in ©-© are 

5B± ■ V ± ~ — — , U± ■ Vj. r — , (A.4) 

£±A A ?±A 

where 9\ is the angle between u± and 5B±, assumed to be small, and • u± = has 
been used to estimate 9\ ~ X/£±\, which is, indeed, small if £j_a ^> A. 

Further development is the same as in the Kolmogorov/GS argument reviewed 
above, except that in (|A.2J) and, consequently, in l)A.3|) . A must be replaced by £j_a- An 
additional assumption is now needed to determine £±\. Boldyrev conjectures that u± 
and 5B± will align to the maximum possible extent. This is achieved if the angle 9\ 
between them is comparable to the characteristic angular wonder of 5B^: 

Aa-t^-T^ => e±A~ (A/||a) V2 . (A.5) 

£-LA ^ || A 

Combining (|A.3jl (where A is replaced by £±\) and (|A.5)I . one gets 

6u±x „ (^) V2 „ (e^) 1 ' 8 „ {£VaX) 1/4 ) (A _ 6) 



The scaling relation (|A.6J) is equivalent to a k 3 ^ 2 spectrum of kinetic energy. 

The status of Boldyrev's theory vis-a-vis real MHD turbulence is uncertain. 
Observationally, only in the SW does one measure the spectra with sufficient accuracy to 
state that they are consistent with k~^ 3 but not with k~ 3 ^ 2 [23121 121] • From numerical 
simulations, it appears that the condition for the k~ 3 ^ 2 spectra [231131 and the alignment 
scaling 6 X ~ (e/v 3 A ) 1/4 X 1 / 4 [H] to emerge is that the mean field is strong (a few times 
5B rms ), 12 whereas in the SW, B ~ SB IIDS . It is not, however, clear why that should 
matter asymptotically, because 8B±\/Bq is arbitrarily small for sufficiently small A. 

12 Note however, that j^Hl had Bq ~ 100SB rms , reported a k ± 3 ^ 2 spectrum, but also found that the 
anisotropy fit the GS scaling lu\ ~ A 2 / 3 , not lu\ ~ A 1 / 2 that appears in l|A.7(l . 



Interplanetary and interstellar plasma turbulence 



15 



Scaling of passive scalar fields 

The scaling of the passively mixed scalar fields, e.g., density fluctuations Sn e , is slaved 
to the scaling of the Alfvenic fluctuations. Again as in Kolmogorov-Obukhov theory, 
one assumes a local-in-scale-space cascade of scalar variance and a constant flux e n 
of this variance. Then, analogously to (jA.ljl . 8n\/r x ~ e n . The cascade time is 
t^ 1 ~ • Vj_ ~ va/1\\ ~ e /& u \\- This gives 

6n x ~ (^y /2 8u ±x , (A.8) 
so the scalar fluctuations have the same scaling as the turbulence that mixes them. 



Kinetic- Alfven-wave turbulence 

The scaling laws for the KAW turbulence are again obtained following the Kolmogorov- 
Obukhov/GS line of reasoning [TJ Locality of interactions and constancy of the 
energy flux imply, analogously to (jA.ljl . 

'^V ^ ~ e B = const. (A.9) 
If the turbulence is strong, then, analogously to (jA.2j) . 



8 5B ± _ 5B X X 

dz B B l\\x 

Assuming that the cascade time is comparable to the inverse KAW frequency, r x ~ 
1\\\X/vaPi, and combining this with (jA.9jl - (jA.10j) . we get 

V3 /„,3\l/3 



SB 



Bn 



A 



—7 /3 

The first of these scaling relations is equivalent to a k ± spectrum of magnetic 
energy, the second quantifies the anisotropy. Note that for KAW-like fluctuations, 
5B\\ X ~ 5B ±X ~ 5B X , 5E X ~ (v APl /c)5B x /X and 5n x /n ~ 55 A /5 (see EZJ. 
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